AD-A162  661  SIMULATIONS  AMONG  MULTIDIMENSIONAL  ITERATIVE  ARRAVS 
ITERATIVE  TREE  AUTONA  (U)  ILLINOIS  UNIV  AT  URBANA 
COORDINATED  SCIENCE  LAB  J  L  TRAHAN  JAN  86 
UNCLASSIFIED  UILU-ENG-86-2282  E/G  12/1 


January  1986 


UILU-ENG-86-2202 

ACT-66 


m 


COORDINATED  SCIENCE  LABORATORY 

College  of  Engineering 


SIMULATIONS  AMONG  MULTI¬ 
DIMENSIONAL  ITERATIVE  ARRAYS, 
ITERATIVE  TREE  AUTOMATA,  AND 
ALTERNATING  TURING  MACHINES 


Jerry  Lee  Trahan 


I 


DTIC 

ELECTEI 
DECS  7  885 


D 


UNIVERSITY  OF  ILLINOIS  AT  URBANA-CHAMPAIGN 


Approved  for  Public  Release.  Distribution  Unlimited. 


85  12  27  02 8 


xil* 


OTITVTV 


!  msnasssnEmEBaam 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

lb.  RESTRICTIVE  MARKINGS 

Unclassified 

None 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  OISTRIBUTION/AV  At  LABILITY  OF  REPORT 

N/A 

Approved  for  public  release,  distribution 

2b.  OECLASSIFlCATION/OOWNG RAGING  SCHEDULE 

unlimited. 

N/A 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER (S) 

6.  MONITORING  ORGANIZATION  REPORT  NUMBER  (SI 

UILU-ENG-86-2202  (ACT-66) 

N/A 

_  6a.  NAME  OP  PERFORMING  ORGANIZATION 

Coordinated  Science  Laboratory 
University  of  Illinois 

6c.  AOORES3  (City.  Slat *  and  ZIP  Coda l 

i’V  1101  W.  Springfield  Avenue 
^  Urbana,  Illinois  61801 

1  Ba.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


Sc.  AOORESS  (City.  Stata  and  ZIP  Coda) 


6.  OFFICE  SYMBOL 
(ItappUe aHa) 

N/A 


8b.  OFFICE  SYMBOL 
(IfappUaaPla) 

N/A 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Office  of  Naval  Research 


7b.  ADDRESS  (City.  Stott  and  ZIP  Coda) 

800  N.  Quincy  Street 
Arlington,  VA  22217 


S.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Contract  #  N00014-85-K-0570 


800  N.  Quincy  Street 
Arlington,  VA  22217 


11.  TITLE  (Inc tad*  Security  Clamtficattoni  5  jjnul a  t  ion S  among 

nal  iterative  arrays,  iterative 


12.  PERSONAL  AUTHOR  IS) 

Trahan,  Jerry  L. 


1 10.  SOURCE  OF  FUNDING  NOS. 

PROGRAM 

PROJECT 

TASK 

WORK  U'-ilT 

ELEMENT  no. 

NO. 

NO. 

NO 

N/A 

N/A 

N/A 

N/A 

13b  TYPE  OF  REPORT 

13b.  TIME  covered 

14.  OATS  OF  REPORT  <Yr..  Mo..  Day) 

is.  page  COUNT 

FROM  TO 

1986  January 

42 

16.  SUPPLEMENTARY  NOTATION 

N/A 

COSATI  COOES 


GROUP  t  SUB.  GR. 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  Mock  number) 

~~We  present  three  simulations:  a  simulation  of  an  alternating  Turing  machine  (ATM) 
operating  in  time  T(n)  by  an  iterative  tree  automation  (ITA)  -4n— time  0(T-frO))>  a  simulation 
of  a  d-d^mensional  iterative  array  (dIA)  operating  in  time  T(n)  by  an  ATM  -in  time- 
0((T(«)-^j£)>,  and  a  simulation  of  an  ITA  operating  in  time  T(n)  by  an  ATM.  in  time’ 

-0( (!!(«■)-)■?}■%  The  first  two  improve  previously  known  results.  The  first  implies  the 
simulation  of  a  nondeterministic  Turing  machine  by  an  ITA  in  time  0(T(n)'>  of  Culik  and 
Yu  (19841i1  The  second  is  stronger  than  the  simulation  of  a  dLA  by  an  ATM  in  tim^ 

0((T(n))  ’  /logT(n))  of  Seiferas  (1977)  and  Dymond  and  Tompa  (1985).  i/  . ,  N 


•* 


20.  0ISTRI3UTION, AVAILABILITY  CF  ABSTRACT 
UNCLASSiFiEO/UNLIMITHO  C  SAME  AS  APT.  Z  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 


22b.  TELEPHONE  NUMBER 
■  Inciuaa  Inro  Codn 


22c.  OFFI 

CS  SYMBOL 

NON 

i _ . 

E 

Unclassified 


MCUHITY  CLASSIFICATION  Of  THIS  PAOi 


11. 


tree  automata,  and  alternating  turing  machines 


SIMULATIONS  AMONG  MULTIDIMENSIONAL  ITERATIVE  ARRAYS. 
ITERATIVE  TREE  AUTOMATA. 

AND  ALTERNATING  TURING  MACHINES 


JERRY  LEE  TRAHAN 


B.S..  Louisiana  State  University  and  A.  &  M.  College.  1983 


THESIS 

Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Electrical  Engineering 
in  the  Graduate  College  of  the 
University  of  Illinois  at  Urbana-Champaign,  1986 

Accession  For 

'  ' "?!  S  GRAM 
r’.T0  T  A  R 
V’  ■K.r-'nrc'.-d 
Ju  ;tivU:jtion_ 


Urbana.  Illinois 


|  Ey - - 

DistrJVntioo/ 

Availability  Codes 
[Avail  and/or 
Dist  I  Special 


Copy  ! 


b&y&'k.M. 


Ill 


abstract 


We  present  three  simulations:  a  simulation  of  an  alternating  Turing  machine  (ATM) 
operating  in  time  Tin  )  by  an  iterative  tree  automaton  (ITA)  in  time  OiTin  )).  a  simulation 
of  a  d-dimensional  iterative  array  (dIA)  operating  in  time  Tin)  by  an  ATM  in  time 
OiiTin)Y).  and  a  simulation  of  an  ITA  operating  in  time  Tin)  by  an  ATM  in  time 
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INTRODUCTION 
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Multidimensional  iterative  arrays,  iterative  tree  automata,  and  alternating  Turing 
machines  are  important  models  of  parallel  computation.  A  d-dimensional  iterative  array 
(dlA)  comprises  an  infinite  set  of  finite  state  machines  located  at  the  points  of  the  d- 
dimensional  integer  lattice.  An  iterative  tree  automaton  (ITA)  consists  of  an  infinite  set  of 
finite  state  machines  connected  into  a  binary  tree.  An  alternating  Turing  machine  (ATM), 
like  a  nondeterministic  Turing  machine  (NTM),  may  have  choices  of  transitions  for  each  com¬ 
bination  of  state,  input  symbol,  and  worktape  symbols.  From  an  existential  stale  an  ATM 
accepts  if  at  least  one  choice  leads  to  an  accepting  state:  from  a  universal  state  an  ATM 
accepts  if  every  choice  leads  to  an  accepting  state.  One  can  view  the  ATM  as  a  computational 
model  that  makes  copies  of  itself  to  evaluate  each  of  the  choices.  A  deterministic  Turing 
machine  (DTM)  has  only  a  single  possible  transition  for  each  combination  of  state,  input 
symbol,  and  worktape  symbols.  Each  of  these  models  has  a  fixed  structure.  Each  processor 
(finite  state  machine  or  ATM  copy)  can  communicate  with  only  a  fixed  set  of  other  proces¬ 
sors.  in  contrast  to  variable  structure  models  such  as  the  hardware  modification  machines  of 
Dymond  and  Cook  (1980)  and  the  parallel  random  access  machines  of  Fortune  and  Wv’lie 
(1978). 

This  thesis  studies  the  simulation  of  a  dlA  by  an  ITA.  The  simulation  of  an  ATM  by 
an  ITA  and  the  simulation  of  a  dlA  by  an  ATM  achieve  this  purpose.  Let  X(t)  denote  the  set 
of  all  languages  recognized  by  machines  of  type  X  in  time  (Ki)  '  here  X  €  {dlA.  ITA,  ATM. 
DTM.  NTMl.  These  simulations  produce  the  following  time  bounds: 

ATM(t)  £  IT  Alt)  and 
dIA(t)  £  ATM  (tJ  ). 

In  addition,  this  thesis  presents  the  simulation  of  an  ITA  by  an  ATM  within  the  bound 
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ITA(l)  £  ATM  (*-).  and  considers  the  simulation  of  an  ATM  by  a  dlA. 

Two  of  the  above  time  bounds  improve  previously  known  results.  The  first.  ATM(t)  £ 
lTA(t).  implies  the  bound  .\TM(t)  £  ITA(t)  established  by  Culik  and  Yu  (1984)  because  an 
XTM  is  more  restricted  than  an  ATM.  The  second.  dIA(t)  £  ATM(t'/  ).  is  better  than  dlA(t) 
£  ATMfr1'  +l/log  t  ).  The  latter  bound  arises  from  the  combination  of  a  result  of  Seiferas 
(1977a).  who  proved  that  dIA(t)  £  DTM(f‘/+1).  and  a  result  of  Dymond  and  Tompa  (1985). 
who  established  that  DTM(t)  £  ATM(t/log  t). 

Aside  from  improving  upon  previous  findings,  the  results  described  in  this  thesis  are 
significant  because  they  extend  current  knowledge  about  ATMs.  In  addition  to  the  usual 
computational  resources  time  and  space  (which  generally  are  inversely  related),  the  ATM  has 
alternations  between  universal  and  existential  states  as  a  computational  resource.  The  ATM 
simulation  of  the  dlA  and  the  ATM  simulation  of  the  ITA  utilize  the  resource  alternation  to 
achieve  the  stated  time  bounds. 

The  ITA  simulation  of  an  ATM  uses  the  ITA's  capability  to  emulate  direct  central  con¬ 
trol.  that  is.  to  act  as  though  the  state  transitions  of  each  finite  state  machine  in  the  ITA 
depend  on  the  state  of  a  uniquely  designated  finite  state  machine.  The  techniques  of  Seiferas 
(1977b)  are  used  to  establish  the  direct  central  control  capability  of  the  ITA. 

The  ATM  simulation  of  the  dlA  uses  the  "divide-and-conquer"  technique,  as  in  Savitch 
( 1970).  Paterson  (1972).  and  Loui  (1981). 

The  framework  of  this  thesis  is  as  follows.  Chapter  2  gives  definitions  of  the  dlA.  the 
ITA.  and  the  ATM.  Chapter  3  reviews  the  literature  related  to  the  thesis  research.  Chapters 
4  and  5  detail  the  simulations  themselves,  prove  their  correctness,  and  determine  their  run¬ 
ning  times.  Chapter  6  describes  the  conclusions  reached  from  this  research  and  offers  some 
open  research  problems. 
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Chapter  2 
DEFINITIONS 

To  more  precisely  describe  the  actions  of  the  computational  models  in  the  simulations 
to  follow,  this  chapter  defines  the  dIA.  the  ITA.  and  the  ATM. 

Let  n  be  the  length  of  the  input  string.  For  all  n,  a  machine  recognizes  an  input  string 
of  length  n  in  time  T(n)  if  the  machine  requires  at  most  T(n)  steps  to  accept  the  input  string 

A  d~dimensional  iterative  array  (dIA)  is  an  infinite  synchronized  d-dimensional  array  of 
finite  state  machines,  one  for  each  d-tuple  of  nonnegative  integers.  In  dIA  M  let  M(X)  denote 
the  machine  corresponding  to  the  d-tuple  X,  called  the  coordinates  of  \1(X ).  Let  0_  denote 
the  d-tuple  of  all  Os.  Only  M(0_)  receives  input  symbols  and  produces  output  symbols.  At 
each  step  M(0_)  reads  a  new  symbol  of  the  input  string  or  does  not  read  a  new  symbol  as  a 
function  of  the  current  state  of  itself,  the  current  states  of  its  neighbors,  and  the  current 
input  symbol.  A  neighbor  of  a  machine  with  coordinates  X  in  the  array  is  a  machine  whose 
coordinates  are  obtained  by  adding  or  subtracting  1  from  a  single  coordinate  of  X.  The  state 
transitions  of  every  machine  depend  on  the  current  states  of  the  machine  itself  and  its  2d 
neighbors.  The  state  transitions  of  M(0_)  depend  additionally  on  the  current  input  symbol. 
A  dIA  accepts  an  input  string  in  I*  when  M(0_)  enters  a  state  in  a  designated  set  of  final 
states. 

Formally,  a  dIA  M  is  a  septuple  M  *  (d  ,Q  X,S((,8,r  ,F )  where 
d  is  a  positive  integer. 

Q  is  a  finite  set  of  states. 

Z  is  a  finite  input  alphabet  ( S  .€£  is  an  endmarker.  B  is  the  blank  symbol ). 

5„:  Q2J  ""‘xi'I  U{5  .$  |  )—Q  is  the  transition  function  for  \1(0_).  the  machine  at  the  origin. 

5:  q2j  *>-.<2  is  the  transition  function  for  every  machine  other  than  M(0  ). 

r.  Q~J  ~'x{Z'dlB  S  1  i—hritc  ,  ‘alse  1  is  the  function  that  specifies  whether  \KO_i  reads  the 
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next  input  symbol,  and 

F  Q  Q  is  a  set  of  final  stales. 

q x  €  Q  is  a  special  quiescent  state.  Initially,  every  machine  is  in  a  quiescent  state.  8 

satisfies  the  condition  8(qk.q  K . q  K)=q K  thus  except  for  M(0.).  a  machine  leaves  the  quiescent 

stale  only  after  a  neighbor  leaves  the  quiescent  state.  Infinite  B's  are  assumed  to  follow  the 
end  of  the  input  string. 

An  iterative  tree  automaton  (ITA)  is  composed  of  synchronized  finite  state  machines 
connected  in  an  infinite  full  binary  tree  structure.  Let  R  (X)  denote  the  root  of  ITA  R.  where 
X  is  the  empty  siring.  In  general,  for  any  finite  binary  string  0  let  R  (00)  and  R  (0 1 )  denote 
the  left  and  right  children  of  R(0).  respectively.  The  Level  of  R(0)  is  1/31.  the  length  of  8. 
Define  I  X  I  =0.  Only  R  (X)  receives  input  symbols  and  produces  output  symbols.  At  each  step 
R  (X)  reads  a  new  symbol  of  the  input  string  or  does  not  read  a  new  symbol  as  a  function  of 
the  current  state  of  the  root,  the  current  states  of  its  two  children,  and  the  current  input 
symbol.  The  current  state  of  the  root,  the  current  states  of  its  two  children,  and  the  current 
input  symbol  determine  the  next  state  of  R(\).  For  every  other  machine  in  the  ITA.  the 
current  states  of  the  machine  itself,  its  parent  and  its  left  and  right  children  determine  the 
next  state  of  the  machine.  An  ITA  accepts  an  input  string  when  R(\)  enters  any  one  of  a 
designated  set  of  final  states. 

Formally,  an  ITA  R  is  a  septuple  R  =  (Q  ,E.8(,.S,  .8,  ,r  .F )  where 
Q  is  a  finite  set  of  stales. 

1  is  a  finite  input  alphabet  ($  €Z  is  an  endmarker.  B  ,€Z  is  the  blank  symbol), 

5,,:  Q  xt  £U  \B  S  !  ixQ:—Q  is  the  transition  function  for  R(K).  the  machine  at  the  root. 

5  .8  C,J— ‘'J  are  the  transition  functions  of  the  left  and  right  children,  respectively .  ol 

each  machine  R  1  3).  3  €  1 0. it*. 

r  <2  x(ZU \B  .S  !  )xQ:  —  {true  .  false  }  is  the  function  that  spec  lies  whether  the  root  reads 
the  next  input  symbol,  and 
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F  £  Q  is  a  set  of  final  states. 

If  X .  W ,  Y .  Z  are  respectively  the  current  states  of  R  (00).  its  parent  £(0).  its  left 
child  £(000),  and  its  right  child  £(001).  then  S,  (X  ,W  .Y  .7.  )  is  the  next  state  of  £(0 0). 
Similarly  8,  specifies  the  next  state  of  each  R  (01). 

€  Q  is  a  special  quiescent  state.  Initially,  every  machine  in  the  ITA  is  in  a  quiescent 
state.  8  satisfies  the  condition  5(q k.q K.q k,q K)=q K  where  8  is  S,  or  8,  thus  except  for  £(A).  a 
machine  leaves  the  quiescent  state  only  after  its  parent  leaves  the  quiescent  state.  Infinite  B's 
are  assumed  to  follow  the  input  string.  The  transition  function  for  a  left  child  must  be 
different  from  the  transition  function  for  a  right  child  for  the  ITA  to  be  distinct  from  a  II A. 

A  configuration  of  an  ITA  R  is  a  pair  (C.w).  where  C  is  a  mapping  from  the  machines  in 
R  to  0  and  w  £L*$  •  For  all  0€{O.1|*.  X  €{0.l}.  where  w  is  the  unread  portion  of  the  input 
string.  R  has  a  Legal  transition  from  (C.w)  to  ( C  \w  ')  if 

(1)  C'(\)=8„(C(\),a.C(0).C(l)). 

(2)  C  ’(0O)=8,  (C  (00 ).C  (000), C  (001  ).C  (0)). 

(3)  C'(01)=8,  (C(01).C(01O).C(011).C(0)).and 

(4)  w  =w  '—a  —$  ,  or 

w  —aw  ’if  r  (C  (\),C  (0),C  ( l),a  )=true  .  or 
w  =w  ’  if  r  (C  (\).C  (0).C  (l).a  )=  false  . 

A  computation  by  ITA  R  is  a  sequence  (C0.w„).  (Cj.wj),  (C;,w:).  ...  of  configurations 
where  the  transition  from  (C;  .wk  )  to  (Ct  ^.^w,  +l)  is  legal  for  all  k.  Note  that  Ct  (0)=C„(0) 
for  ;  0  '  fi-k  since  we  stare  with  all  machines  in  a  quiescent  state. 

An  alternating  Turing  machine  (ATM)  is  defined  in  Chandra,  kozen,  and  Stockmeyer 
(1981).  A  configuration  of  an  ATM  is  the  state  of  the  ATM.  the  contents  of  the  input  tape, 
the  contents  of  each  of  the  worktapes.  and  the  locations  of  each  of  the  tape  heads.  An  ATM 
is  a  Turing  machine  in  which  every  nonfinal  state  is  either  universal  or  existential  A 
configuration  with  an  existential  state  is  accepting  if  at  least  one  successor  configuration  is 
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accepting.  A  configuration  with  a  universal  state  is  accepting  if  every  successor  configuration 
is  accepting.  An  ATM  accepts  an  input  string  if  its  initial  configuration  is  accepting.  An 
ATM  has  a  two-way  read-only  input  tape  with  endmarkers  and  k  worktapes.  which  are  ini¬ 
tially  blank  A  step  of  an  ATM  consists  of  reading  one  symbol  from  each  worktape  and 
reading  an  input  symbol,  then  writing  a  symbol  on  each  of  the  worktapes.  moving  each  of 
the  heads  left  or  right  one  tape  square  or  not  moving  the  tape  heads,  and  choosing  a  new  state 
from  the  set  specified  by  the  transition  function. 

One  can  describe  all  possible  computations  of  an  ATM  on  some  input  siring  as  a  compu¬ 
tation  tree.  All  nodes  are  configurations,  the  root  is  the  initial  configuration,  and  the  children 
of  any  configuration  c  are  exactly  those  configurations  that  can  be  reached  from  c  in  one  step 
according  to  the  transition  rules  of  the  ATM.  The  leaves  of  the  tree  are  the  final 
configurations  and  may  be  accepting  or  rejecting.  A  branch  of  the  computation  tree  is  a 
downward  directed  path  from  the  root:  in  other  words,  a  branch  is  a  sequence  of 
configurations  starting  with  the  initial  configuration.  Assume  that  for  an  ATM  to  run  in  time 
T(n).  all  branches  terminate  in  at  most  T(n)  steps. 

Formally,  an  ATM  AM  is  a  septuple  AM  -  ( k  ,Q  .Z.T.8.q„.g  )  where 

k  is  the  number  of  wotkcapes. 

Q  is  a  finite  set  of  stales. 

£  is  a  finite  input  alphabet  (S  /€£  is  an  endmarker). 

T  is  a  finite  worktape  alphabet  (B  €  T  is  the  blank  symbol). 

S:  Q  XT  *  X(£  U|5  ;)— * ‘P(Q  X(T  —{B  1)‘  x{left  ,  right .  stationary  i4  *')  is  the  transition 
function,  where  P(S)  is  he  power  set  of  S.  that  is.  the  collection  of  subsets  of  S. 

q„  k  Q  is  the  initial  state,  and 

g:  Q  —{universal  .  existential ,  accept  .  reject  )  is  a  mapping  identifying  each  stale  as  a 
universal,  existential,  accepting,  or  rejecting  state. 
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Chapter  3 

LITERATURE  REVIEW 


This  chapter  reviews  the  literature  related  to  the  research  reported  in  this  thesis. 

Chandra.  Kozen.  and  Stockmeyer  (1981)  present  the  concept  of  alternation.  (The  same 
authors  originally  presented  the  concept  in  Chandra  and  Stockmeyer  (1976)  and  Kozen 
(1976).)  This  thesis  uses  their  ATM  model.  They  derive  significant  relationships  between 
classes  of  languages  accepted  by  time  and  space  bounded  ATMs  and  those  accepted  by  time 
and  space  bounded  DTMs.  In  particular,  logarithmic  alternating  space  is  equivalent  to  poly¬ 
nomial  deterministic  time,  and  polynomial  alternating  time  is  equivalent  to  polynomial 
deterministic  space.  Paul.  Prauss.  and  Reischuk  (1980)  demonstrate  that  an  ATM  with  a  sin¬ 
gle  tape  can  simulate  an  ATM  with  multiple  tapes  in  linear  time. 

Dymond  and  Tompa  (1985)  prove  another  result  related  to  this  thesis.  They  establish 
that  DTM(t)  £  ATM(t/log  t).  Their  proof  associates  the  computation  of  the  DTM  with  an 
acyclic  directed  graph.  They  use  a  two-person  pebbling  game  to  pebble  the  graph  within  a 
time  bound  of  (Xn/log  n)  for  a  graph  with  n  vertices.  Next,  the  ATM  steps  simulate  the 
two-person  pebbling  of  the  graph.  In  the  pebbling  game,  one  person's  moves  correspond  to 
existential  choices  of  the  ATM.  and  the  other  person's  moves  correspond  to  universal  choices 
of  the  ATM. 

Paterson  (1972)  represents  a  TM  computation  as  a  two-dimensional  diagram  of  succes¬ 
sive  tape  configurations.  He  employs  divide-and-conquer  in  both  time  and  space  dimensions. 
This  method  is  generalized  in  this  thesis  in  the  simulation  of  a  dlA  by  an  ATM.  Loui  (1981) 
establishes  a  space  bound  for  a  DTM  to  accept  the  same  language  as  a  d-dimensional  \TM 
with  one  worktape  head.  The  proof  utilizes  a  generalization  of  crossing  sequences  across  the 
boundaries  of  d-dimensional  boxes  of  the  worktape.  He  uses  a  divide-and-conquer  method  to 
recursively  partition  the  boxes. 
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Rosenfeld  (1979)  presents  a  good  review  of  iterative  automata. 

Cole  (1969)  formally  presents  the  d-dimensional  iterative  array  of  finite  state 
machines.  He  establishes  that  the  computing  speed  of  a  dIA  can  be  increased  by  a  constant 
factor  by  enlarging  the  set  Of  states  of  each  machine.  He  proves  that  the  class  of  context-free 
languages  does  not  contain  all  the  sets  of  strings  accepted  by  a  dIA.  nor  do  the  sets  of  strings 
accepted  by  a  dIA  contain  all  context-free  languages.  He  proves  that  computing  capability 
increases  as  the  number  of  dimensions  increases. 

Seif  eras  (1977a)  extends  Cole's  work  on  deterministic  dIAs  to  nondeterministic  dlAs 
(NdlAs).  He  derives  that 

NTMU' )  Q  NdIA(t). 

NdIA(t)  C  NTM(r'  *0.  and 
dIA(t)  C  DTM(frf*0. 

The  second  result  is  related  to  the  simulation  of  an  NdIA  by  an  ATM  given  in  this  thesis. 
His  simulation  uses  about  nd  steps  of  a  one-dimensional  (2d+l)  head  TM  to  simulate  the  nth 
step  in  a  computation  of  a  dIA. 

Seiferas  (1977b)  establishes  that  a  dIA  with  direct  central  control  is  no  more  powerful 
than  a  regular  dIA.  and  that  a  regular  dIA  can  simulate  a  dIA  with  direct  central  control  in 
linear  time.  In  this  simulation,  the  finite  state  machine  at  the  origin  of  the  dIA  controls  the 
dIA  indirectly  by  propagating  the  value  of  its  state  outward  using  only  local  communication. 

Culik  and  Yu  (1984)  construct  a  language  L  such  that  an  ITA  accepts  L  in  real-time, 
but  no  dIA  can  accept  L  in  real-time.  They  state  that  the  converse  problem,  that  is,  whether 
real-time  dIA  languages  ( d  ^2)  are  properly  contained  in  real-time  ITA  languages,  is  an  open 
problem.  They  establish  that  an  ITA  can  simulate  an  NTM  in  linear  time.  Their  simulation 
provides  a  basis  for  the  simulation  of  an  ATM  by  an  ITA  in  Chapter  4. 


Chapter  4 

THE  ITA  SIMULATION  OF  THE  DIA 
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This  chapter  contains  a  simulation  of  the  ATM  by  the  ITA  and  a  simulation  of  the  dIA 
by  the  ATM  and  proofs  of  the  correctness  of  each  simulation. 

Lemma  1:  Every  language  recognized  in  time  t  by  a  k  -tape  ATM  can  be  recognized  in  time 
O  (t )  by  a  one- tape  ATM. 

This  result  is  from  Paul.  Prauss.  and  Reischuk.  (1980).  They  specify  that  the  one-tape 
ATM  does  not  have  separate  input  and  output  tapes. 

Lemma  2:  Every  t  steps  of  an  ATM  with  at  most  c  ^  3  choices  at  each  step  can  be  simulated 
by  (c  —1  )t  steps  of  an  ATM  with  at  most  2  choices  at  each  step. 

Proof:  Let  M  -  ( k  ,Q  .1  .r.8jj0.g  )  be  an  ATM  with  at  most  c  choices  at  each  step.  We 
define  an  ATM  M '  =  ( k  .Q  '.L.r.8'.?0  .  g)  with  at  most  2  choices  at  each  step  such  that  M' 
simulates  M.  Let  a  transition  rule  of  M  be 

S^o.w^a,,)  =  {(</,.w1.X,),(y2.»v2.X2) . ( qj  .Wj  ,X;  )| 

for  q{).q  t,  .  .  .  ,q1  €Q  .  w  0.w  , . w  €T*  ,  a,)6ZUj5  I.  Xi,X2 . Xy  €\left  .right  .stationary  \. 

O^/'^c.  If  y'^2.  then  8'(9o.w0.a0)sS((y0,iy0,a0).  Otherwise.  S'  has  the  following 

corresponding  rules: 

S‘{q  n.w  ,,.anj  -  iq  ,.w  ,.X  ,).(p  x.w stationary  )} 

S'(p  |.W(,.a„)  =  ;(y2.w2,X2)  (p2.w„.jfaftnnarv  i| 

5 '(p  _2.w,,.<2,,;  =  {(</._l.iv._l.Xy  _,).(</,  ,w ,  X  !. 

The  new  state  set  Q‘  includes  Q  and  all  the  new  slates  p..p2 . p  _2  that  are  used  in  the 
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decomposition  of  the  steps  with  three  or  more  choices  into  steps  with  two  choices. 

It  is  clear  that  each  move  of  M  can  be  simulated  by  M'  with  at  most  c— 1  moves.  So 
(c  —1  )r  steps  of  M '  are  enough  to  simulate  t  steps  of  M.  □ 

Definition:  Subtree  Broadcast  ITA  (SBITA)  —  An  SBITA  S  is  the  same  as  an  ITA  except 
as  specified  in  the  following.  Some  of  the  machines  in  S  are  designated  to  be  control  units.  If 
S(a)  is  a  control  unit,  then  for  every  descendant  S(0)  of  S(a).  the  next  state  of  S(0) 
depends  on  the  current  state  of  S(a)  as  well  as  on  the  current  states  of  S(0)  and  its  parent 
and  its  children.  At  all  times,  for  every  machine  S  (0).  at  most  one  ancestor  of  S  (0)  is  a  con¬ 
trol  unit.  During  the  computation  a  control  unit  S(a )  ceases  to  function  as  a  control  unit 
when  it  designates  its  children  S(aO)  and  S(erl)  to  become  control  units.  Initially.  S(A)  is 
the  only  control  unit. 

Lemma  3:  Every  t  steps  of  an  SBITA  can  be  simulated  by  It  steps  of  an  ITA. 

Proof:  Let  S  be  an  SBITA.  We  design  an  ITA  R  that  simulates  S. 

R  (0)  is  designated  as  a  control  unit  when  the  state  of  R  (0)  is  in  a  special  set  of  states 
called  the  broadcast  set.  The  broadcast  set  contains  a  particular  state  qpass  .  and  R  (0  )  enters 
state  qpa„  on  the  step  before  it  will  pass  its  control  unit  capabilities  to  its  descendants.  If 
R(/S)  is  in  state  qp„ss  at  any  time  t.  then  at  time  t+1.  R(&)  will  be  in  some  state  not  in  the 
broadcast  set.  At  any  time  t.  the  stale  of  R  (0)  can  be  in  the  broadcast  set  only  if  one  of  the 
following  is  true: 

( 1 )  0=A  and  t-<). 

(2)  at  time  t-1.  the  state  of  #(0)  was  in  the  broadcast  set.  or 

(3 )  at  time  t- 1 .  the  parent  of  R  ( 0 )  was  in  state  . 

The  remainder  of  the  proof  is  taken  directly  from  Seiferas  (1977b).  and  is  modified  to 
apply  to  deterministic  IT  As  rather  than  nondetermimstic  dIAs.  Informally,  a  machine  in  R 
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simulating  a  control  unit  in  S  performs  a  broadcast  by  propagating  the  values  of  its  states  to 
its  descendants.  The  simulation  of  each  descendant  will  lag  by  time  equal  to  its  distance 
from  the  control  unit,  so  communication  between  a  control  unit  and  any  descendant  will 
take  twice  as  long  in  R  as  in  S. 

Let  Q  '  denote  the  set  of  states  of  each  machine  in  S.  Let  S'0.  8‘, ,  and  S',  denote  the  tran¬ 
sition  functions  of  S. 

Let  Q  denote  the  set  of  states  in  R.  where  q  \=({0.1)x^x^x«/;J. 
Q  =(\OA)xQ'XQ'xQ‘)Uq A.  When  q  €Q  .  denote  the  first  component  as  phase(q).  the  second 
component  as  control(q)  the  third  component  as  prev(q),  and  the  fourth  component  as 
current(q).  Control  (q  )  propagates  the  states  of  the  control  unit  to  its  descendants:  previq  ) 
holds  the  previous  state  of  a  machine  for  reference  by  its  children:  current  (q  )  holds  the 
current  state  in  the  simulation  of  the  corresponding  machine  in  S. 

Define  for  R  (A)  as  control  unit 

(0,  8 '0(control  (qx  ). current  {qx  ).a  .q .q  ).  current  (y,  ). 

8'Q(control  (qx  ), current  (qx  ).a  .q  ,q\  )) 
if  €Q-<i  \-phase(qx  )*1. 

( 1  .control  (qx  ),prev  (qx  ). current  (qx  )) 

8„(^v  .a  .q,  .q,  )  =  if  qx  SQ-q^.  phase  (qx  )=0.  r  iqx  ,a  ,qt  .q,  )=  false  . 

(1  .phase  (qx  ), phase  (<?.,  ),phase{qx  )) 
if  V,  =?.v  r  (</,  -a  .qt  .q,  )=  false  ,  or 

0 

otherwise  . 


phase  (q,)  if  q,  sq  \  or 
where  q  ,  —  Cltrrgnl  )  otherwise. 

Define  for  R  (0).  where  R  (0)  is  a  control  unit. 


12 


5/  (y.  4t  4,  4P  )  = 


(0.  S'/  (.control  (q,  ). current  (qy  ).q  jq  ‘r  jq  'p  ),  current  (qy  ). 

S',  (control  (qx  ).current  (qx  ).q .q  .q  'p  )) 
if  qy  £Q—qvphase(qx  )=1. 

(1. control  (qx  ).prev(qx  ). current  (qx  )) 
if  qx  £Q  —q  \.phase  (qx  )=0. 

( 1, phase  (qx  ). phase  (qx  ).phase(qx  )) 

>f  ?.«•  =9  or 


0 

otherwise  . 


where 


phase  (q,  )  if  q,  =q  x  or 
current  (qt)  otherwise. 


S,  is  denned  similarly. 

Define  for  R  (0).05*\.  where  R  (0)  is  noi  a  control  unit. 


8/  (?.v  4t  4,  4P  )  = 


(0. control  (qp  ) .current  (qx  ).  S',  (control  (qx  ).current  (qx  ).q  .q  .q  'p  )) 
if  </v  4?  €<2  ~<i  s.phase  (qx  )=1 . 

(Note :  This  transition  propagates  control  information.  ) 

( 1  .control  (q x  ),prev(qx  ).currenl  (qx  )) 
if  qA.phase(qx  )=0. 

( 1  .control  (qp  ). phase  (qx  f. phase  (qx  )) 
if  ?>  -q.\4P  \.or 

0 


otherwise. . 


where 


h\- 


phase  (q,  )  if  q,  =q  \ 

prev  (q,  )  if  q.  €Q  —q  vi  —p  .  or 

current  (q,)  otherwise. 


5,  is  defined  similarly. 


Informally,  the  first  three  cases  in  each  definition  above  are  "simulate  a  transition." 
"wait  your  turn  to  simulate  a  transition,"  and  "set  up  to  start  simulation."  respectively. 


This  completes  the  definition  of  IT  A  R  that  simulates  SB1TA  S. 
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A  configuration  of  R  is  a  pair  (C  .w  );  a  configuration  of  S  is  a  pair  (C  \w  ). 

The  sequence  a=[(C  „.w  ,(C  t.w  ,  ).(C2.^2)..--1  of  configurations  of  R  is  skewed  if  ail  of 
the  following  hold  for  all  nonnegative  integers  j  .  0  €  (0.1)*: 

w2j  =w2/+l- 

if  j  <  10  I .  then  Cj  (0)=C(l(0)=y  v 
if  j  >  101.  then  C;  (0)€(2  . 

1  for  y  — 1  /3  I  <’dd  , 
phase  (C;(0))=  0  for  y._,0, 

even  . 

control  (C  131+;  +1(0))=currenr  (C;  +1(\)). 
prev(C  ,0i+|(/3))  -  current  (C  i<jh-i(/3)) 

»  prev(C  t3i+2(0)) 

-  /wev  (C  isi+jW) 

=  p/iare  (C, >(0)).  and 

current  (C  131+2;  +3(3))“  currenl  (C  131+2/  +3W) 

-  prev(C  i3i+2;  +«(0)) 

-  previC  101+2/ +j(3))- 

For  such  a  skewed  sequence  a.  define 

skew -1(a)=((C  \},w  '(>).(C'i.w  'jUCS.w  '2 )....] 
if  the  following  hold  for  all  nonnegative  integers  j  .  0  €  {0,1  (*.  q  -C  131+2/  (0): 

w  ' .  =w  2;  and 

phase  (q)  dq—q\or 
currenl  (q  )  otherwise 

The  definition  of  a  skewed  sequence  does  not  restrict  the  choice  of  values 
current  (C{  4,-2, (0))  for  j>  0.  so  every  sequence  (C’,.lvV-  (C  Vw of 

configurations  of  S  where  all  machines  are  initially  quiescent  is  equal  to  for  some 

skewed  sequence  a.  We  show  below  that  a  is  a  computation  by  R  from  (C„.w„'i  if  and  onlv 
if  skew~Ha)  is  a  computation  by  S  from  (C„.w,J.  It  will  suffice  to  map  state  q  6(2  to  state 


C',(3)= 
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I 

phase  (q  )€Q  'if  q—q  A  or  to  state  current  )€Q  ’  otherwise.  j 

i 

Let  a  be  a  skewed  sequence  with 

< 

a=[(C0.w()).(C  i.v*’0).(C2.w  j).(C3.w  and  | 

skew  ~‘(a)=[(C  '0.w0).(C  'i  .w  j).(C  '2.w  2)....].  , 

To  prove  that  a  is  a  computation  by  R  if  and  only  if  skew~l{at )  is  a  computation  by  S.  it 

suffices  to  verify  for  all  nonnegative  integers  j,  0  €  10.1  }*  U  {X}.  y  =  I  j8 1  or  y  *  1 0  I  +2  j  +2. 

that  C>v+|(0)  results  from  the  transition  from  l.C(j3,  y)  and  that  C  t^i+2;+2(0)  results  from 

the  transition  from  LC(f3. 1  0  I  +2j  +1 )  if  and  only  if  C'y+1(0)  results  from  the  transition 

from  LC(jS.j).  □ 

Lemma  4:  Every  t  steps  of  a  one-head,  one-tape  TM  can  be  simulated  by  a  two-stack 
machine  in  up  to  5t  steps. 

Proof:  This  lemma  is  proved  in  Hopcroft  and  Ullman  (1979). 

Call  the  stacks  of  the  two-stack  machine  stack  A  and  stack  B.  Let  tp(l) .  tp(j) 

denote  the  nonblank  contents  of  the  TM  tape.  The  two-stack  machine  represents  a  TM 

configuration  by  hoiding  tp(l) . tp(i)  in  stack  A  with  tp(i)  at  the  top  and  tp(i+l) . tp(j) 

in  stack  P  with  tp(i+l)  at  the  top.  where  tp(i)  is  the  tape  element  that  the  head  of  the  TM  is 
reading. 

In  five  steps,  the  two-stack  machine  can  simulate  a  TM  step  in  which  the  head  moves  to 
the  righ 

1  ’op  tp(i)  from  stack  A. 

2.  Change  the  symbol  contained  in  tp(i)  to  the  symbol  to  be  written  in  this  step. 

3.  Push  tp(i)  back  onto  stack  A, 

4.  Pop  tp(  i-1 )  from  stack  B. 

5.  Push  tpt  i-rl )  onto  stack  A. 


In  three  steps,  the  two-stack  machine  can  simulate  a  T\1  step  in  which  the  head  moves 
to  the  left: 

1.  Pop  tp(i)  from  stack  A. 

2.  Change  the  symbol  contained  in  tp(i)  to  the  symbol  to  be  written  in  this  TM  step. 

3.  Push  tp(i)  onto  stack  B. 

In  three  steps,  the  two-stack  machine  can  simulate  a  TM  step  in  which  the  head  remains 
on  the  same  tape  element: 

1.  Pop  tp(i)  from  stack  A. 

2.  Change  the  symbol  contained  in  tpCi)  to  the  symbol  to  be  written  in  this  TM  step. 

3.  Push  tp(i)  back  onto  stack  A.  □ 

Theorem  5:  For  all  T(n  ).  every  language  recognized  in  time  T(n  )  by  a  k  -tape  ATM  can  be 
recognized  by  an  ITA  in  time  O  ( T  ( n  )). 

Proof:  Let  AM  be  a  one-tape  ATM  with  at  most  two  choices  at  each  step.  By  Lemmas  1, 
2.  and  3.  to  prove  the  theorem  it  suffices  to  show  that  an  SBITA  S  can  simulate  AM  in  time 
0(T(n)). 

It  is  clear  that  a  one-head,  one-tape  TM  can  directly  implement  a  branch  of  the  compu¬ 
tation  tree  of  AM.  By  Lemma  4,  a  two-stack  machine  can  simulate  the  actions  of  a  one-head, 
one-tape  TM. 

S  can  implement  the  stacks  in  the  manner  described  below  (modeled  after  the  imple¬ 
mentation  in  Culik  and  Yu  (19K4)).  For  all  |3  and  all  x>0.  all  descendants  of  -S'(j3)  at  dis¬ 
tance  a  from  S(0)  are  in  the  same  state  and  contain  the  same  stack  elements.  Let  A(3) 
denote  the  portion  of  S  consisting  of  S(i 3)  and  all  of  its  descendants.  A(/3)  implements  a  pair 
of  stacks,  stack  A  and  stack  B.  each  of  which  operates  separately  according  to  the  rules 
below.  Every  machine  in  A (3)  is  a  finite  state  machine  and  can  store  one  stack  element  of 
each  stack.  S(t 3t  is  a  control  unit  and  holds  the  state  of  the  two-stack  machine  and  the 
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elements  at  the  top  of  each  stack.  To  push  stack  A.  S(0)  broadcasts  to  each  of  its  descen¬ 
dants  5(a)  to  send  the  element  of  stack  A  contained  in  S' (a)  to  both  children  of  S' (a),  and 
5 (0)  stores  the  element  to  be  pushed.  To  pop  stack  A.  S (0)  broadcasts  to  each  of  its  descen¬ 
dants  5  (a)  to  send  the  element  of  stack  contained  in  S  (a)  to  the  parent  of  S’  (a). 


At  the  beginning  of  the  simulation,  by  definition  of  an  SBITA.  S(A)  is  the  only  control 
unit.  A(A)  is  functioning  as  a  pair  of  stacks.  A  and  B.  to  store  the  input  string  w.  S(A)  reads 
w  and  stores  w  in  stack  A.  5(A)  next  pops  w  from  stack  A  and  pushes  w  onto  stack  B.  At 
this  point,  w  is  stored  in  stack  B  with  the  first  symbol  of  w  at  the  top  of  stack  B.  In  addi¬ 
tion.  A(A)  is  functioning  as  a  pair  of  stacks.  C  and  D.  to  simulate  the  actions  in  a  branch  of 
the  computation  tree  of  AM.  This  simulation  continues  as  long  as  each  step  in  the  computa¬ 
tion  tree  of  AM  has  only  one  choice.  When  AM  first  takes  a  step  with  two  choices,  control 
unit  5(A)  broadcasts  to  its  descendants  to  push  all  stacks  one  level  down,  and  S(A)  passes 
the  control  unit  capability  to  5(0)  and  5(1).  Then  A(0)  simulates  one  of  the  choices,  and 
All)  simulates  the  other  choice.  5(A)  no  longer  acts  as  a  control  unit,  but  enters  existential 
state  </,  if  AM  is  in  an  existential  state  at  this  point  or  enters  universal  state  qu  if  AM  is  in  a 
universal  state  at  this  point.  Whenever  two  choices  are  present  at  a  node  of  the  binary  com¬ 
putation  tree.  S  will  proceed  to  simulate  each  of  the  choices  in  the  manner  described  above 
except  the  steps  will  be  for  A(0)  instead  of  for  A(A).  S  simulates  each  branch  of  the  binary 
computation  tree  in  parallel  with  the  other  branches. 

When  a  branch  finishes  a  computation,  the  unique  control  unit  for  that  branch  sends  the 
answer  (acceptance/ rejection )  up  to  its  parent.  5(0).  0€’().1|*.  If  5(0)  is  in  state  jr  .  it  will 
enter  state  y,,ic  .  indicating  acceptance,  or  enter  stale  q,r.  .  indicating  rejection,  based  on  the 
table  in  Figure  4.1;  if  5(0)  is  in  state  u,.  .  it  will  enter  state  <v„.  or  slate  </, . ,  based  on  the 
table  in  Figure  4.2. 

To  determine  the  time  required  for  S  to  simulate  AM.  look  at  the  various  components 
of  the  simulation.  The  reading  of  input  string  w  and  positioning  of  w  in  slack  B  in  5(A) 


\  -v-\w 


requires  time  2n.  By  Lemma  4.  the  two-stack  simulation  of  each  branch  requires  time 
0(T(n)).  The  concurrent  SB1TA  simulation  of  two-stack  machines  requires  time  T(n).  The 
passing  of  the  result  up  to  the  root  requires  time  T(n).  Therefore  the  overall  simulation  of 
AM  by  S  requires  time  0(T(n)).  □ 


Next  is  the  ATM  simulation  of  the  dIA. 

Theorem  6:  For  all  T(n  ).  every  language  recognized  in  time  Tin  )  by  a  dIA  can  be  recognized 
in  timeO((7’(n  ))‘/ )  by  an  ATM. 

Proof:  Let  M  be  a  nondeterm inistic  dIA  of  time  complexity  T(n).  We  design  an  ATM 
AM  that  will  simulate  the  operation  of  M.  (Note:  AM  will  existentially  guess  T(n).) 

The  local  configuration  of  M(2L)  at  lime  t.  denoted  LC(X.t).  is  the  (2d+l)-tuple  contain¬ 
ing  the  stale  of  M(X)  at  time  t  and  the  states  of  the  2d  neighbors-  of  MCXj.)  at  time  t. 

The  computation  of  M  on  the  input  can  be  expressed  by  a  (d+1  )-dimensional  computa¬ 
tion  array  C  in  which  the  value  at  coordinates  (X,t)  is  the  state  of  M(X)  at  time  t.  except  for 
one  column  of  C  that  contains  the  input  string.  (See  Figure  4.3.)  For  all  d-dimensional 
subarrays  D  of  C  denote  the  endpoints  of  D  as  x  Io  and  x  lb  in  dimension  1 .  x2„  and  c  2b  in 

dimension  2 .  xJt,  and  xJh  in  dimension  d,  and  ta  and  th  in  dimension  d+x.  where 

x  \a  <*i*.  x2a  <x2 b . la  <*h  ■  For  k=l d,  let  denote  U*,,—  ljcu, Xu>  ,xtb  +1 1.  Let 
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Let  QD  be  the  boundary  of  D.  $D  contains  a  value  for  each  coordinate  vector  in  the  uni 
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Informally.  QD  consists  of  6 d  +2  sets  of  values  around  the  boundary  of  D.  The  values  at  the 
corners  of  D  are  contained  in  multiple  sets.  Each  set  of  values  is  contained  on  a  separate  tape. 
QC  contains  the  input  siring. 

(Note:  For  d-1.  QD  is  shown  in  Figure  4.4.) 

Assume  that  when  MfO^)  enters  an  accepting  state,  then  on  the  following  steps  V1(0_) 
remains  in  the  accepting  state,  and  every  other  machine  in  M  proceeds  to  enter  the  quiescent 
state. 

Define  valid  (QD  )  to  be  a  predicate  that  is  true  if  and  only  if  there  is  an  assignment  of 
states  to  the  entries  of  D  such  that  all  state  transitions  of  each  machine  in  M  follow  from  the 
transition  rules  of  M.  Call  such  an  assignment  consistent .  The  recursive  procedure  ASIMD. 
defined  below,  computes  the  predicate  valid  .  that  is.  ASIMD  returns  "true"  or  "false"  ASIMD 
uses  a  divide-and-conquer  method  according  to  lime  and  according  to  each  dimension  of  the 
machines  in  M. 


Initially.  AM  reads  the  input  string,  existentially  inserts  null  symbols  in  the  input 
string  for  the  steps  when  M  does  not  read  an  input  symbol,  existentially  guesses  the 


machines 


time 
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t* 


Figure  4.4  -  QD  ,  input  to  procedure  ASIMD  for  1-dimensional  case 


remainder  of  QC ,  and  calls  AS1MD(§C ).  AM  assumes  that  the  length  of  every  side  of  C  is 
the  smallest  power  of  2  greater  than  or  equal  to  T(n).  AM  accepts  the  input  string  if  and 
only  if 

( 1 )  ASIMD(QC  )  returns  "true."  and 

(2)  the  stale  of  M(0.)  at  time  T(n)  determined  by  ASIMD($C  )  is  an  accepting  state. 

The  recursive  procedure  ASIMIX^D  )  is  defined  as  follows: 

Input:  QD  . 

Input  invariant:  There  exists  an  integer  j  such  that  for  all  i  x,h—xia  -2J  and  th—ta  =2J 
Call  2J  the  width  of  D. 

Output:  "true"  or  "false" 

Step  1:  If  the  width  of  D  is  1.  then  check  that  for  all  Y  representing  coordinates  in  D. 
the  state  of  M(Y)  at  time  th  =r„  +1  correctly  follows  from  LC(Y.f„  )  according  to  the  transi¬ 
tion  rules  of  M.  If  this  condition  holds,  then  the  answer  returned  by  this  invocation  of 
AS1MD  is  "true";  otherwise,  the  answer  is  "false." 

Step  2:  At  this  point  the  width  of  D  is  greater  than  1.  Let  t,„  =(f„  +th  )/2.  Existentially 

guess  the  states  for  all  machines  whose  coordinates  are  in  ( X\ . X'd )  at  time  t,„ .  For  all 

machines  whose  stale  at  time  t,„  is  in  QD  .  check  that  the  state  guessed  in  this  step  is  equal  to 
the  state  in  QD  .  If  the  states  are  not  equal,  then  return  "false." 

Step  3:  For  each  i.  1  .  let  x„„  =(x,u  +x.h  )/2.  Existentially  guess  LC(Y.t)  for  all  Y 

€  (.Y  i  ,X2 . .t„„ . Xd  )  and  all  t„  < i  ,. .  Check  that  the  states  in  l.C(Y_  .t ,  ).  l.C  ( F  .1 ..  ). 

and  LC(Y.t,„)  are  ;qual  to  the  corresponding  stales  in  QD  and  to  the  corresponding  states 
guessed  in  Step  2.  If  any  of  the  corresponding  stales  are  not  equal,  then  return  ’’false."  (See 
Figure  4.5.) 

Step  4:  Universally  choose  one  of  the  2J  ~l  subarrays  of  D  created  from  performing 
div  ide-and-conquer  in  each  of  the  (d+1)  dimensions  oi"  D. 
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Step  5:  Let  S  be  the  subarray  chosen  in  Step  4.  Call  AS1MD(Q5  ). 

Step  6:  If  all  universal  choices  in  Step  5  return  "true."  then  return  "true";  otherwise, 
return  "false." 

Now  we  prove  that  AM  accepts  if  and  only  if  M  accepts. 

First  we  show  that  ASIMD  (QC  )  returns  "true"  if  and  only  if  valid  (QC  )  is  true.  Let  2J 
be  the  width  of  QD  .  the  input  to  ASIMD.  By  induction  on  j  we  show  that  ASIMD  (QD  ) 
returns  "true"  if  and  only  if  valid  (QD  )  is  true. 

Basis  (j-1):  Step  1  of  ASIMD  confirms  that  each  state  at  time  th  follows  from  the 
states  at  ttl  according  to  the  transition  rules  of  M.  Steps  2  and  3  confirm  that  state  informa¬ 
tion  on  separate  tapes  for  the  same  machine  at  the  same  time  is  equal.  Therefore,  valid  ($D  ) 
is  true. 

Conversely,  if  valid  ($D  )  is  true,  then  the  assignment  of  states  to  entries  of  D  Q$D  is 
consistent.  Therefore.  Step  1  of  ASIMD  confirms  that  each  state  at  time  th  follows  from  the 
states  at  ta  according  to  the  transition  rules  of  M.  and  Steps  2  and  3  confirm  that  state  infor¬ 
mation  on  separate  tapes  for  the  same  machine  at  the  same  time  is  equal.  Therefore. 
ASIMD  ($D  )  returns  "true." 

Inductive  hypothesis:  ASIMD  (QD  )  returns  "true"  if  and  only  if  valid  ($D  )  is  true. 

Induction  (j>l):  Assume  ASIMD  (QD  )  returns  "true."  Let  E.  F . G  be  the  subarrays 

of  D  from  Step  4  of  ASIMD.  It  follows  that  the  calls  to  ASIMD  {QE).  ASIMD(QF) . 

ASIMD  (QG  )  returned  "true."  By  the  inductive  hypothesis,  valid  {^E).  valid  (QF) . 

valid  (Q/G  )  are  all  true.  Since  array  D  is  completely  covered  by  the  overlapping  subarrays  E. 

F . G.  and  the  stales  assigned  to  overlapping  locations  of  E.  F . G  are  the  same  since  they 

are  passed  to  ASIMD  (QE  ).  ASIMD  ($F  ) . ASIMD  (QG  >  from  ASIMDiQD  ).  and  there  is  a 

consistent  assignment  of  states  to  the  entries  of  E.  F . G.  then  there  is  a  consistent  assign¬ 

ment  of  states  to  the  entries  of  D.  Therefore,  valid  (0D  '  is  true- 
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Conversely,  assume  valid  (QD)  is  true.  It  follows  that  valid  (QE).  valid  (QF) . 

valid  (QC  )  are  all  true.  By  the  inductive  hypothesis.  ASIMD  (QE).  ASIMD  (QE) . 

AS/MD  (0G  )  all  return  “true."  Therefore,  by  Step  6  of  AS1MD.  ASIMD  (QD  )  returns  "true." 

Now  we  show  that  if  AM  accepts,  then  M  accepts.  If  AM  accepts  input  string  w.  then 
ASIMD  (3C  )  returns  "true.“  and  the  state  of  M(0  )  at  time  T(n)  determined  by  ASIMD  ($C ) 
is  an  accepting  state.  Since  ASIMD  ($C )  returns  "true."  valid  (&C )  is  true.  Because 
valid  (0C )  is  true  and  the  state  of  M(0)  at  time  T(n)  determined  by  ASIMD  (0C  )  is  an 
accepting  state,  an  assignment  of  states  to  all  machines  in  M  for  all  t  .  O^t  <T(n  ).  that  ends 
with  M(0.)  in  an  accepting  slate  exists,  such  that  the  states  of  all  machines  at  time  t  result 
from  the  states  of  all  machines  at  time  t  —1  according  to  the  transition  rules  of  M  and  input 
string  w.  This  implies  that  M  accepts  w. 

Next  we  show  that  if  M  accepts,  then  AM  accepts.  If  M  accepts  input  string  w.  then 
there  exists  an  assignment  of  states  to  all  machines  in  M  for  all  t.  0 ^T(n).  such  that  the 
states  of  all  machines  in  M  at  time  c  follow  from  the  states  of  ail  machines  at  time  t  -1 
according  to  the  transition  rules  of  M  and  the  input  string  w.  and  such  that  all  machines  are 
initially  quiescent,  and  such  that  M(0 )  is  in  an  accepting  state  at  time  T(n).  As  a  result. 
valid  (qC  )  is  true:  hence.  ASIMD  (QC  )  returns  "true."  Because  ASIMD  (QC  )  returns  "true." 
and  M(0.)  is  in  an  accepting  state  at  time  T(n).  AM  accepts  w. 

Now  we  show  that  the  time  required  for  this  simulation  is  0{{T in  )Y  ).  Let  D  denote 
a  ( d+ 1  )-diraensional  subarray  of  length  k  in  each  dimension.  Let  TM,{k  )  denote  the  lime 
complexity  of  the  simulation.  ASIMD  (QD  )  is  performed  on  computation  array  D  which  has 
sides  of  length  k  .  Then  AM  selects  one  of  the  subarrays  S  of  D  and  calls  ASIMD  (^S  ).  The 
lengths  of  the  sides  of  S  are  k  /2.  so  ASIMIXQS  )  requires  time  Tnlik  /2).  The  lime  to  per¬ 
form  Steps  2  and  3  of  ASIMD  {QD  )  is  0  ( k  J  ).  The  time  complexity  of  ASIMD  ($D  )  is 

TM,(k  )  =  TAV(k  !2)+OikJ  )  =  0<kJ  ). 

The  space  required  for  this  simulation  is  Oik1'  )  In  particular,  the  simulation  of  M  by  AM 


takes  time  7\*A/  (T(n))  -  0((T(n  )Y  )  and  space  O  (( Tin  )Y  ).  □ 


The  same  techniques  can  be  used  to  simulate  an  NdlA  on  an  ATM  in  the  same  time. 

Corollary:  For  all  T(n  ).  every  language  recognized  in  time  7*(n  )  by  an  NdlA  can  be 
recognized  in  time  O  (( T (n  )Y  )  by  an  ATM. 


Chapter  5 

THE  ATM  SIMULATION  OF  THE  IT  A 


This  chapter  contains  a  simulation  of  the  ITA  by  the  ATM  and  a  proof  of  the  correct¬ 
ness  of  the  simulation.  It  also  outlines  a  simulation  of  an  ATM  by  a  dlA. 

Theorem  7:  For  all  Tin  ).  every  language  recognized  in  time  Tin  )  by  an  ITA  can  be  recog¬ 
nized  in  time  OiiTin  ))2)  by  an  ATM. 

Proof:  Let  R  be  an  ITA  of  time  complexity  T(n).  We  design  an  ATM  AM  with  four 
worktapes  that  will  simulate  the  operation  of  R.  (Note:  AM  will  existentially  guess  T(n).) 

For  0€{O.1|*  and  X  €<0,l).  the  local  configuration  of  Ri&X)  at  lime  l.  denoted 
LC  (0X  x  ).  is  the  quadruple  iqy  (t  (t  ).q,  it  ).qp  it  ))  where 
q x  it  )  is  the  state  of  R  (0X  )  at  time  t. 
q i  i t )  is  the  state  of  R  (/3X  0)  at  time  t. 
q,  (t )  is  the  state  of  R  it 3X  1 )  at  time  t. 

qp  it )  is  the  state  of  R  (0)  at  time  t  if  0X  or  qp  (r  )  is  the  input  symbol  at  time  t  if 
0X=A. 

Define  I  \  I  =0. 

Let  A (0)  denote  the  portion  of  R  that  comprises  Ri0)  and  all  descendants  of  Ri 0).  Let 
<r, 3  denote  the  sequence  of  states  qy  it  )  of  R  (0)  for  O^t  ^Tin  ).  Define  valid  (cr 3V  cr3,0X  ) 
to  be  a  predicate  that  is  true  if  and  only  if  there  is  an  assignment  of  states  to  A (0X  )  for 
every  lime  t  ,  O^t  ^  Tin  ).  such  that  for  every  i  .  0^ i  ^  T in  ).  the  stales  of  every  machine  in 
A(0X  )  at  lime  i  follows  from  the  states  of  every  machine  in  A(0X  )  at  time  t— 1  according 
to  the  transition  rules  of  R.  crav .  and  cr3.  Call  such  an  assignment  of  states  consistent.  The 
recursive  procedure  ASIMT.  defined  beiow.  computes  the  predicate  valid  .  that  is.  ASIMT 
returns  " true "  or  " false" 
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Initially.  AM  writes  input  siring  w  onto  tape  4.  existentially  inserting  null  symbols  for 
the  steps  when  M  does  not  read  an  input  symbol.  AM  then  existentially  guesses  the  states 
for  £(X)  at  each  lime  step  from  0  to  T(n).  then  calls  ASIMT (<rK.w  A).  AM  accepts  if  and 
only  if 

( 1 )  the  output  of  ASIMT (o\.w  A)  is  "true."  and 

(2)  the  state  assigned  to  R  (X)  at  time  T(n)  by  ASIMT  (o\.w  A)  is  an  accepting  state. 

The  recursive  procedure  ASIMT (<Tp v  ,(Tp.0X  )  is  defined  below. 

Inputs:  c Tpx  (the  sequence  of  states  qy  (t  ))  on  tape  1.  <7S  (the  sequence  of  states  qn(t  )) 
on  tape  4.  and  0X . 

Output:  "true"  or  "false" 

Step  I:  If  I  0X  I  —Tin  ).  then  return  "true"  if  all  stales  of  </>  it  )  are  quiescent:  if  not  all 
states  are  quiescent,  then  return  "false."  If  !  0X  I  <Tin  ).  then  perform  the  following. 
Existentially  guess  <rPXl,  on  tape  2.  Existentially  guess  crp v,  on  tape  3.  (See  Figure  5.1.) 
Verify  for  each  t  .  0<t  ^ Tin  ).  that  qv  (t )  follows  from  LC{04  —1).  defined  on  tapes  1-4. 
according  to  the  transition  rules  of  R.  If  any  </x  U  )  does  not  follow,  then  return  "false."  Ver¬ 
ify  that  the  symbols  generated  for  qp  it )  are  equal  to  the  symbols  already  on  tape  4.  If  they 
are  not.  then  return  "false."  If  no  consistent  sequence  of  guesses  is  possible,  then  return 
"false." 

Step  2:  Copy  the  contents  of  tape  1  onto  tape  4. 

Step  3:  Universally  choose  either  R  (0X  0)  or  Rit 3X  1).  If  R(0XO)  is  chosen,  then  copy 
the  contents  of  tape  2  onto  tape  1.  then  call  ASIMT (<rax ».<Tax  .0X0).  If  Ril 3X  1)  is  chosen, 
then  copy  the  contents  of  tape  3  onto  tape  1 .  then  call  ASIMT  (<7a.v  \.crax  ,0X  1 ). 

Step  4:  If  both  universal  choices  in  Step  3  return  "true."  then  return  "true":  otherwise 
return  "false." 
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Now  it  will  be  shown  that  AM  accepts  if  and  only  if  R  accepts. 

We  show  by  induction  on  h  —Tin. )— !  j3X  I .  that  ASIMT  (crux  .cr3./3X  )  returns  "true"  if 
and  only  if  valid  (o^y  ,<r^.j3X  )  is  true. 

Basis  (h-0):  If  ASIMT (crgV  .crfl./3X  ).  I  0X  I  =T  in  ).  returns  "true."  then  the  sequence  of 
local  configurations  is  consistent  with  all  descendants  of  Ri j3X )  according  to  the  transition 
rules  of  R.  since  all  machines  in  A(/3X  )  are  in  a  quiescent  state  at  all  times  and  have  no  effect 
on  the  computation.  This  implies  that  valid  (cr3.y  ,cr3.|3X  )  is  true. 

Conversely,  if  valid  (<r3>  .(T^.fiX  )  is  true,  then  a  consistent  assignment  of  states  to 
A(0X  )  exists.  According  to  the  transition  rules  of  R.  all  states  assigned  to  A10X  )  must  be 
quiescent,  because  all  machines  in  A(/3X  )  are  below  level  T(n)  in  R.  Therefore.  Step  1  of 
ASIMT  id 3S  .(Tfj.fiX)  confirms  that  all  states  in  cr3y  are  quiescent.  Therefore. 
ASIMT  icrax  .cr3.0X  )  returns  "true." 

Inductive  hypothesis:  ASIMT  id Sx  .cra.3X )  returns  "true"  if  and  only  if 

valid  (cr3V  .cra.3X  )  is  true. 

Induction  !h>0):  Assume  ASIMT  (o^y  0’a-&X  )  returns  "true."  It  follows  that  for  some 
crayo  and  (Tpxi  the  calls  to  ASIMT  (<t  ax  H’CTqx  &X  0)  and  ASIMT  (<r3.y  [.cr^y  ,/3X  1)  returned 
"true."  By  the  inductive  hypothesis,  valid  (o^.yo.o^.y  .j3X  0)  and  valid  (<r3y  j.Ujy  .&X  1 )  are 
both  true.  By  Step  1  of  ASIMT id3x  ,c Tj.QX  ).  cr3  v  is  a  sequence  of  states  q,  it  )  for  R(3X  ). 
such  that  (/,  it  )  follows  from  LC(<3X ./  —1)  according  to  the  transition  rules  of  R.  Therefore, 
a  consistent  assignment  of  states  to  A(3X  )  exists.  Therefore,  vrlid  (c r3V  ,cr3.3X  )  is  true. 

Conversely.  assume  valid  <.o"3y  ,cr3,3X )  is  true.  It  follows  that  both 
valid  (<t3x <i&a.\  3X  0)  and  valid  (o,fl.y  ,,<r3-V  .0X  1 )  are  true.  By  the  inductive  hypothesis, 
both  ASIMT ( cr3y  ,,.<T3y  ,3X  0)  and  ASIMT  (<T3.y  ]  ,<Tfl y  ,j3X  l )  return  ’’true.”  Therefore,  by  Step 
4  of  ASIMT.  ASIMT (<t3x  ■O'e-QX  )  returns  "true." 

\ow  we  show  that  if  AM  accepts,  then  R  accepts.  If  AM  accepts  input  string  w.  tnen 
ASIMT((TK,w  A i  has  returned  "true."  and  the  stale  assigned  to  R(K)  at  time  T(n  '  by 


ASIMT  (a^.w  A)  is  an  accepting  state.  Because  ASIMT i(JK.w  A)  returns  "true," 
valid  ( o\,w  A)  is  true.  Since  valid  (o\.w  A)  is  true,  and  since  the  state  assigned  to  R  (A)  at 
time  T(n)  is  an  accepting  state,  R  accepts  w. 

Next  we  show  that  if  R  accepts,  then  AM  accepts.  If  R  accepts  input  string  w.  then 
there  exists  an  assignment  of  states  to  all  machines  in  R  for  all  i  .  O^t  (n  ),  such  that  the 
states  of  all  machines  in  R  at  time  t  follow  from  the  states  of  all  machines  at  time  r  —1 
according  to  the  transition  rules  of  R  and  input  string  w.  and  such  that  all  machines  are  ini¬ 
tially  quiescent,  and  such  that  R(\)  is  in  an  accepting  state  at  time  T(n).  As  a  result. 
valid  (crv.w  A)  is  true.  Since  valid  ( o\,w  A)  is  true.  ASIMT  {crK.w  A)  returns  "true."  Because 
ASIMT (cr\.w  A)  returns  "true,"  and  the  state  of  R  (A.)  at  time  T(n)  is  an  accepting  state.  AM 
accepts  w. 

Now  we  show  that  the  time  required  for  this  simulation  is  0((T(n  ))2).  Let  TAAl(h  ) 
denote  the  lime  complexity  of  ASIMT  (c rfl  v  .ap.&X  ).  where  h  stT  in  )—  I  j3X  I .  Steps  2  and  3 
of  ASIMT (cr^v  .ar^fiX  )  require  time  O(Tin)).  ASIMT (cr j3_v  ,cr 0./3X  )  calls 
ASIMT  {<Tqxy  -Vqx  .&XY  ),  Y  €  {0.1 1,  which  has  a  time  complexity  of  TA\t  (h  —  1 ). 
HO)  =  OiTin  )).  These  terms  give  rise  to  the  recurrence 

Tm/CH  )=TAA/(h  ~l)+0(T(n  )). 

In  particular,  the  simulation  of  T(n)  steps  of  R  by  AM  requires  time 

TAlW(T(n))=TA,w(T(n)-l)+0(T(n))=0((T(n))2). 

This  simulation  requires  space  0  (T  ( n  )).  □ 

In  order  to  comolete  the  cycle  of  simulations  from  IT  A  to  dIA  and  from  dIA  to  IT  A. 
only  the  simulation  of  an  ATM  by  a  dIA  remains.  In  order  to  obtain  a  time  bound  on  this 
simulation,  first  consider  a  simulation  of  an  NTM  by  a  dIA.  Suppose  that  a  (deterministic) 
dIA  operating  in  time  O  ).  where  p  is  some  constant,  could  simulate  i  i  steps  of  an  NTM. 
Seiferas  (1077a)  proves  that  a  DTM  operating  :n  time  0(i<  *'  !  can  simulate  : ;  steps  of  a 
dIA.  Together,  these  two  simulations  would  imply  that  a  DTM  could  simulate  an  NTM  in 
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polynomial  time;  hence,  the  computational  complexity  class  P  would  equal  the  computational 
complexity  class  NP.  The  equality  P  =  \P  is  widely  believed  to  be  unlikely;  hence,a  polyno¬ 
mial  time  simulation  of  an  NTM  by  a  dlA  is  unlikely.  Therefore,  a  polynomial  time  simula¬ 
tion  of  an  ATM  by  a  dlA  is  unlikely. 

Proposition  8:  For  all  T(n).  every  language  recognized  in  time  T(n)  by  an  ATM  can  be  recog¬ 
nized  in  time  0{ 2ri"  *)  by  a  dlA. 

Briefly,  a  (deterministic)  dlA  can  simulate  an  ATM  in  exponential  time  as  follows.  The 
dlA  can  simply  compute  each  branch  of  the  compulation  tree  of  the  ATM  in  turn.  Hach 
branch  corresponds  to  the  actions  of  a  DTM.  and  a  dlA  is  able  to  simulate  a  DTM  in  linear 
time  according  to  Seiferas  ( 1977b). 

Theorem  7  and  Proposition  8  together  yield  an  exponential  time  bound  for  a  dlA  simu¬ 
lation  of  an  1TA.  One  would  expect  this  bound  because  the  number  of  finite  state  machines 
potentially  involved  in  a  compulation  grows  polynomiallv  with  time  for  a  dlA.  but  grows 
exponentially  with  time  for  an  ITA. 
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Chapter  6 

CONCLUSIONS  AND  OPEN  PROBLEMS 


This  thesis  has  presented  three  simulations  and  discussed  a  fourth.  When  combined. 
Theorems  5  and  6  imply  that  an  ITA  can  simulate  a  dIA  in  time  OUJ  ),  and  Theorem  7  and 
Proposition  8  imply  that  a  dIA  can  simulate  an  ITA  in  exponential  time. 

These  results  and  the  work  done  in  obtaining  them  suggest  several  open  problems. 

1.  Can  the  time  bounds  of  Theorems  5.  6,  and  7  be  improved? 

2.  Is  there  a  language  L  such  that  some  dIA  recognizes  L  in  linear  time,  but  every  ITA 
requires  superlinear  lime  to  recognize  L?  Culik  and  Yu  (1984)  pose  this  question,  but  for 
real-time.  One  candidate  considered  for  L  is  a  string  of  the  form 

/,  =(*,#*  ,#...#*,„  ##yI#y2#...#v„,  ). 

where  .t  ,..t2 . .r„,  is  an  unordered  list  of  values,  and  y  j.y2 . y,„  is  a  sorted  list  of  the  same 

values.  This  candidate  fails  because  an  ITA  can  sort  in  lime  0(n)  according  to  Browning 
(1979).  and.  though  a  2IA  can  sort  in  time  Oi'/n  logn  )  according  to  Thompson  and  Rung 
(1977).  Nassimi  and  Sahni  (1979).  and  Stout  (1982).  the  dIA  must  write  down  the  output, 
leading  to  a  total  time  requirement  of  0(n).  A  second  possible  candidate  is  a  language  of 
binary  strings  that  represent  connected  d-dimensional  figures. 

3.  How  much  time  is  required  for  an  X-tree  array  to  simulate  an  ATM?  An  X-tree  is  a 
binary  tree  with  additional  edges  connecting  all  nodes  at  the  same  level  in  the  tree.  An  X- 
tree  array  is  an  iterative  array  of  finite  state  machines  organized  into  an  X-tree. 

4.  How  much  time  is  required  for  an  ATM  to  simulate  an  X-tree  array? 

5.  How  can  an  ITA  with  depth  as  a  function  of  the  length  of  the  input  siring  simulate 


an  ATM  or  an  NTM? 
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6.  How  much  time  and  space  are  required  for  an  ATM  with  a  limit  on  the  number  of 
its  alternations  to  simulate  either  the  dlA  or  the  ITA? 

A  further  possible  area  for  future  research  is  the  alternating  iterative  array,  such  as 
either  an  alternating  dlA  or  an  alternating  ITA.  The  addition  of  universal  choices  to  non- 
deterministic  dIAs  or  ITAs  adds  a  second  kind  of  parallelism.  Initial  work  could  be  done  in 
relating  such  a  model  to  other  models  of  computation. 
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